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Abstract 
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constraints from positivity and unitarity. 
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1 Introduction 



In many respects, the two-Higgs-doublet model (2HDM) provides the simplest extension of the Stan- 
dard Model (SM) Higgs sector. For a general introduction to the model the reader is referred to [1], or 
to the more recent papers [2-5] dealing with various aspects of the 2HDM. The extensive review [6], 
discussing the role of the 2HDM in the MSSM, might also be useful. 

Despite its simplicity, the general 2HDM is an interesting theoretical laboratory, adding possible 
new phenomena such as spontaneous and explicit CP-violation, lepton number violation, charged 
Higgs bosons, a dark matter candidate, and a heavy neutral Higgs boson compatible with electroweak 
precision data. Since the 2HDM suffers the same quadratic divergences in corrections to the scalar 
masses as the SM, these precise theoretical constructs of the 2HDM are perhaps not likely to be 
manifested in a more complete theory of nature - and certainly not all of them at the same time. 
However, the large theoretical interest in these models calls for the development of phenomenological 
tools to help theorists provide experimentalists with viable models. This is the purpose of the two- 
Higgs-doublet model calculator (2HDMC) which we describe here. 

2HDMC (Two-Higgs-Doublet Model Calculator) is a C++ code which can be used to study the phe- 
nomenology of a general 2HDM. At the moment, only the CP-conserving 2HDM has been implemented. 
The Higgs potential in 2HDMC can be specified in different parametrizations, allowing for everyone to 
use their own favorite. The Yukawa sector can be specified with full freedom - including specifying 
couplings which lead to fiavour-changing neutral currents. Some standard choices (like type I, type H 
etc.) can also be used. The code provides methods to check theoretical properties of the model, such 
as positivity of the potential and 5-matrix unitarity. There are also methods for calculating the 2HDM 
contribution to the oblique electroweak parameters {S,T,U,V,W,X), to determine the contribution 
to {g — 2)^, and to check the model against constraints from colliders. In addition, decay widths of 
the Higgs bosons can be calculated. The whole code is fully modular and based on object-oriented 
principles. It can be used as a library, in stand-alone applications, or exchanging data with other 
codes through a LesHouches-inspired format. 

The organization of this paper is as follows: section [2] contains a short, hopefully self-contained, 
overview of the physics of the 2HDM. The focus here lies, of course, on the aspects relevant for 2HDMC. 
Section [3] discusses theoretical and experimental ways of constraining the 2HDM parameter space, 
while section [3] describes the structure of the 2HDMC code itself and how to use it. How to obtain and 
compile the code is described in section [5l If you are already familiar with the 2HDM and want to 
get started right away, you should probably skip straight to section HI 

2 Physics of the 2HDM 
2.1 The Higgs potential 

Denote a complex scalar S'C/(2)l doublet field Assigning a hypercharge Y = 1, the component 
representation of <I> becomes <I>(x) = {^'^ (x) , {x))^ . In its most general form, the 2HDM contains 
two identical such fields <I>i and related by a global U{2) symmetry under which $a — > Uab^b and 
$i ^\ul^. In the conventions of [4], the most general gauge invariant and renormalizable potential 
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that can be formed using two Higgs doublets is given by 



Vgen =m^i$|$i + 771^2^2^2 " mi2^i^'2 + h-C 



+ ^Ai + ^A2 {^l^2y + A3 (^l^-l) ($^2^2) + A4 ($1«'2) f^-^^l 

+ Ua5 {^\<^2y + [A6 + A7 {^1^2)] («'I^2) + h.C 



The parameters m^j^, 77x22 ^^'^ Ai_4 are real numbers (since the potential has to be real), whereas 
the remaining parameters A5_7 and 771^2 ™ general can be complex. Non-zero imaginary parts of 
the complex parameters which cannot be removed by a rephasing transformation give rise to explicit 
CP-violation in the Higgs sector. Since we are not going to treat CP-violating effects here, we assume 
all parameters to be real in the following. 

Since the potential given by eq. ([T]) is manifestly [/(2)-invariant, values specified for the parameters 
{mfj,Xi} can only have definite (physical) meaning when a particular basis is specified for the scalar 
fields. Alternatively, one may reformulate the 2HDM using a fully basis-invariant language [4] (see 
also [7] for a group theoretic analysis of the basis invariance). Since this formalism has so far been 
sparingly applied in 2HDM phenomenology, we will maintain the notion of choosing a particular basis. 
Expressed in terms of vacuum expectation values (vevs), a generic basis respecting the {7(1)em gauge 
symmetry can be written as 



_v_ / 

^ V cos/3 



i^2) = — ( / 

V2 V e*^sin/3 



where v = {V^Gi 



-1/2 



246 GeV. By convention 0</3<7r/2is chosen. A non-zero phase ^ results 



in a vacuum which breaks CV spontaneously. Since we do not aim to discuss CP-violating effects, we 
take ^ = 0. Note that the value of tan/3 = ($2) / (^1) at this point only determines one particular 
choice of basis. Since the 2HDM potential is invariant under a change of this basis, tan (3 can not be 
a physical parameter of the model in general. Within the set of CP-conserving bases defined above 
there exists a special choice, called the Higgs basis, in which only one of the two doublets is assigned 
a non-zero vev. 



2.2 Electroweak symmetry breaking 

The electroweak SU{2) x C/(l) symmetry manifest in eq. ([T]) is spontaneously broken down to U{1)em 
by a negative eigenvalue of the scalar mass matrix mfj , causing at least one of the Higgs doublets to 
develop a vev. In a basis where < tan/? < 00, the two minimization conditions 

mfi = 777^2 tan /3 — -v'^ ^Ai cos^ /? + (A3 -|- A4 + A5) sin^ (3 + SAe sin /? cos P + Xr sin^ (3 tan (3^ , (2) 

771I2 = "^12 cot P — -v"^ 1^X2 sin^ /3 -|- (A3 + A4 -|- A5) cos^ P + Xq cos^ P cot P + 3A7 sin P cos /3^ (3) 

can be used to eliminate and 777,22 from eq. Similar equations for the minimum conditions, 
not given here, exist for the Higgs basis. Eliminating and 777-22 leaves eight real parameters in the 
Higgs potential, not counting tan/3 and v. 
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When the symmetry is broken, the eight degrees of freedom from the SU{2) doublets $i and $2 
are usuahy re-expressed in states with definite physical properties. The spectrum then contains three 
Goldstone modes: and G'^ , which are absorbed to give mass to the gauge bosons and Z, 
reducing the number of physical Higgs states to five. Three of these states are neutral, of which two 
{h and H, with rrih < mn) are CP-even, and one (typically denoted A) is CP-odd. The remaining 
two are a pair of charged Higgs bosons (H^). The charged Higgs bosons are a characteristic of the 
2HDM not present in the SM. 

Introducing the mixing angle a to diagonalize the mass matrix for the CP-even states, the original 
doublets are expressed as 

^, = J-( ^/2(G+cos/3-F+sin/3) \ 

y/2 V ^ cos (5 — h sin a + H cos a + i (G'^ cos (3 — A sin /3) J 

^ \/2 (G+sin/3 + F+cos/3) \ 

The Higgs masses can now easily be calculated for any choice of Higgs potential and basis. Using the 
generic potential of eq. ([1]), one arrives at the relations [2] 



2 2 

= . "T^^ a - ^ (2-^5 + As cot p + Xj tan/?) 
sm p cos p 2 



,2 



(6) 



mjj+ = y (A5 - A4) . 

Since the CP-even states mix, their mass matrix M is given by 



where 



jg2 _ f Aic| + 2\QSpcp + X^sj^ (A3 + A4) spcp + Aec^, + Ays^ > 
~ V (A3 + A4) spcp + A6c| + Ar^l A2s| + ^X-jSpcp + X^c?^ J ' 

Performing a rotation by a to diagonalize M, one obtains the masses of the CP-even states 

""o^ ^2)=na)M'n^{a) (9) 

with > m|. Explicitly, the mass eigenvalues are given by [2] 

1 " 



Mi, +Mh±J {Mj, - Ml,f + 4 {Ml,)' 



(10) 



Because of gauge invariance, and the way masses has to be supplied to the EW gauge bosons, 
all couplings of Higgs bosons to gauge bosons are completely determined in terms of the invariants 
Sf3-a = sin(/3 — a) and c^_q = cos(/3 — a). For this important angle we find it convenient to use a 
convention in which —tt/2 < (3 — a < it/2. This choice allows the basis to be changed by varying 
tan/3 while preserving the invariants S/j— a and c/3-a- It is important to respect this convention when 
working explicitly with the angle a, e.g. when interfacing 2HDMC with other codes. 
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2.3 The Yukawa sector 



In the following we will assume the same fermion content as in the Standard Model (with massless 
neutrinos). The mass eigenstates for the down- and up-type quarks are written as vectors in flavour 
space, denoted D and U respectively, and similarly for the leptons L and z^. The 2HDM Yukawa sector, 
expressed in the physical Higgs mass eigenstates {h, H, A, H^), is then given in basis- independent form 
as [4] 

-^Yukawa =^=d\^K^ Sp^^ + p'^ Cp^^^Dh + -^D^Hi^Cp-a - p^Sp.a]DH + -^D-f^p^DA 

+ l=u{K^sp.a + p^cp.^]uh + 1=u[k^cp_^ - p^sp_^]uH - ^Uj^P^UA 

+ p{VcKMp''PR - p^VckmPl}DH+ + up'^PnLH+ + h.c. 

The diagonal 3x3 matrices are given by = \f2M^ jv^ where are the corresponding mass 
matrices for the fermions [F = D,U,L), and as usual Pji/l = (1 i 75)/2. The generality of this 
Yukawa Lagrangian is present through the freedom to choose the components of each p^ arbitrarily. 
A symmetric p^ , as required by CV conservation, has six independent components. In practice, 
the allowed size of the off-diagonal elements in p^ are strongly constrained, since non-zero elements 
induce Higgs-mediated flavour-changing neutral currents (FCNC) at tree-level. Dealing with this 
problem through a careful adjustment of individual elements in p^ involves severe fine-tuning. From 
a theoretical standpoint, this can therefore be deemed an unnatural solution. 

To address this problem in a more elegant way, one can invoke the theorem of Glashow and 
Weinberg [8], stating that tree- level FCNC are absent in any theory where a given fermion does not 
couple to more than one Higgs doublet. This situation can be arranged by imposing a suitable Z2 
symmetry to distinguish $1 from <I>2- Assigning definite Z2 quantum numbers also to the right-handed 
fermions, the undesired couplings can be avoided. When all fermions of a given type (up-quarks, 
down-quarks and leptons) have one common Z2 quantum number, the combinatorial possibilities for 
the Yukawa couplings are enumerated in 2IIDM language as types I, II, III, and IV [9]|l| If the Z2 
symmetry is to be explicit in the Yukawa sector, a dependence is introduced on the basis chosen for 
the Higgs sector. This way tan f3 is promoted to a physical parameter through its appearance in the 
Yukawa couplings. In table [U we list the relations between p^ and for the physically distinct types 
of Yukawa couplings realized with a Z2 symmetry. We note finally that some of the terms present in 
the Higgs potential violate this symmetry. The Ag and A7 terms in particular lead to hard Z2 violation 
and should therefore be suppressed, while an term in many cases can be tolerated since it only 
violates Z2 softly (meaning that the symmetry is restored in the UV). This issue is discussed in some 
detail in [3]. 

2.4 Decay widths 

From the masses and couplings, decay widths of the Higgs bosons {Hi = h, H2 = H, H-^ = A, 
H4 = H^) into various modes are calculated. At the tree-level, the Higgs bosons may decay into pairs 

*Type III is often used in different context referring to a more general class of models with a broken (e.g. by SUSY 
loop effects) Z2 symmetry in the Yukawa sector. In such models tan /3 is not a physical parameter. 
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Type 






I 


II III 


IV 




cot (3 


— K^tan/? — K^tan/J 


cot /3 




cot [3 


cot P cot /3 


cot /? 




cot /? 


— K^tan/3 cot P 


— tan P 



Table 1: Relations between Yukawa coupling matrices and in the four different types of 2IIDM 
Yukawa sectors where all fermions of a given type {F = U, D, L) share the same Z2 quantum number. 



of fermions, pairs of gauge bosons, one gauge boson and another Higgs boson, or into pairs of lighter 
Higgs bosons. The decays Hi — > 77 and Hi ^ gg (z = 1 . . . 3) are induced at the one-loop level. 

2.4.1 Hi ^ //' 

We denote the coupling of Higgs boson Hi [i = 1 .. .A) to a pair of fermions by Cs for the scalar 
part, and by Cp for the pseudoscalar part. As an example, we get from eq. (jlip that the couplings 
of hdd are Cs = [K^sp-a + p^cp-a) and Cp = 0. Only for the charged Higgs are both C5 and 
Cp (potentially) non-zero at the same time. Introducing kinematic variables xi = mf/mHi, and 
X2 = nip/ ruH^ , the decay width to fermions is given by 

r(F, ^ //') = [1 - {x, + x,f] + [1 - (xi - x^f] |CpP}aV2 (1, , (12) 

where as usual 

K'^,x,y) = {l-x-yf -Axy. (13) 

In eq. (|12p . Nc refers to the number of "colour" degrees of freedom for the given fermion, i.e. Nc = 2> 
for quarks and A'^c = 1 for leptons. For the decay of a neutral Higgs into quarks, the QCD radiative 
corrections are known at order [10]. When these are included, the corrected width reads 



1 + 5.67^ii^ + (35.94 - 1.36n,) ( 



(14) 



in the MS scheme with Uf light flavours. The Ug correction is applied also to the decay of the charged 
Higgs. For both neutral and charged Higgs decays, resummation of the leading logarithmic corrections 
is implemented to all orders by using the running MS fermion masses mf{p) in the Higgs couplings. 
The renormalization scale is then set at the Higgs mass scale fi = rriHi ■ 



2.4.2 Hi VV 

With k = roy /rnj^, , the on-shell decay width of a neutral Higgs boson i^j (« = 1 . . . 3) to a pair of 
massive gauge bosons is given by 

IzSvr my 

The fact that two Z bosons in the final state are identical is handled by the factor by = 2 (1) for 
V = W {Z). Note that the coupling CniW in our convention has the dimension of mass, which makes 
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the overall dimension of this expression check out. When m//. < 2my, the three-body decay with one 
off-shell gauge boson decaying into a pair of fermions can instead be considered. Following [11], we 
express the Dalitz density in the scaled fermion energies xi^2 = 2-E'i,2/"T'Hi as 

dr . (1 - xi)(l - X2) + ^(2x1 + 2x2 - 3 + 2k) 

[Hi^vv V jij2) = ^HiVv vm ' 



d3;idx2 ' (1 - xi - X2)2 + k^v 

where the scaled width is = ry/ml^.. Integrating between the kinematic boundaries 

k 

1 — X2 — k<Xi<l 



1 - X2 (17) 
< X2 < 1 - /c 



the off-shell width is obtained. The normalization factor is given by 



Kh.vv = -^^^\CH,vv\'^rnH,S'y, (18) 



with = 1 for the W coupling to fermions, and 



6'z = ^-^sm'9w + ^sm'ew. (19) 



2.4.3 Hi VHj 



It is also possible to have the Higgs Hi (i = 1 ... 4) decay to a combination of one Higgs- (j = 1 . . . 4) 
and one massive gauge boson. The width in this case becomes 



ICnVHl'^rnh ( "^h" "^if- \ 1/0/ 

nn, - VH,) = mS^A 1, -f , -f^ A^/^ 1, 



m. 



IGvr^ niHi \ ' fny ' my j \ ' m?^ ' m?^ 



(20) 



where V denotes one of the gauge bosons Z, , or W~ , such that the decays Hi — > W~H~^ and 
Hi W~^H~ are distinct. In addition to the on-shell decay, it is interesting to consider the decay 
when the gauge boson is off-shell. Analogous to eq. (jl6p . the Dalitz density in the rescaled energies 
xi and X2 of the final state fermions can be expressed [12] as 

jjy, Hf\L)=KHVH - - X2) - kH .21) 

dxidx2 * ^ ^ {I - xi - X2 - kn + ky Y + ky^v ' 

Here kn = iT^\./m\., ky = my/m'jj., and jy = Ty/m'jj^. The integration boundaries are given by 



eq. (jl7p . with k replaced by ky. The overall normalization factor is 



Khvh = ^^^\^H,yHfml6'y (22) 



where 5[y = 1 and 6'^ is given by eq. (fT9l) . 
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2.4.4 Hi HjHk 

For the decay to a pair of scalars (i = 1 ... 3, j,k = 1 ... 4) we only consider the on-shell width, 
since the three-body decay width depends on the subsequent decay modes of the Higgs bosons. With 
Xj = m'jj./m'jj., and Xk = 'fn'^H /m^j., the width is given by 



Note that 5jk = is used for the decay into a pair of charged Higgs bosons. 



(23) 



2.4.5 Ha 
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The two-photon decay of the neutral Higgs bosons Hi [i = 1 ... 3) proceeds through a loop with 
contributions from fermions, W bosons, and charged Higgs bosons. The width is known [13,14], and 
can be summarized in the expression 



V{Hi ^ 77) 



■Ml 



2567r3i 



^(|5^(MhJ|2 + |P^(MhJ|2), 



(24) 



where the form factors are 



S^{Mh,) = 2j2NcQ}C^„jf-—Fs{Tf) - Ch,ww7^F,{tw) - CH,H+H-7r^Fo{rH+) (25) 



for the scalar amplitude, and 



f 



ruf 



(26) 



for the pseudoscalar part. As usual, Nc = 3(1) for quarks (leptons), and Qf is the electric charge in 
units of the elementary charge |e|. The various functions appearing in these form factors are related 
in terms of a common scaling function /(r): 



F,(T)=T-i[l + (l-r-i)/(T) 
Fo{t)=t-' [r-V(T)-l] 



Fpir) = 



Fi{t) = 2 + 3t-^ + 3t-\2 - T-^)f{T) 



(27) 



The argument is r/ = mjj./Amj for the fermions and equivalently for the and H^ contributions. 
Finally, the (complex) function /(r) itself is given by the integral 



fir) = -l r^ln[l-4ry(l-y)] = <( 

^ Jo y 



arcsin^ (v^) 



In 



Vt + Vt-1 



-1 2 



ITT 



T < 1 



T > 1. 



(28) 



2.4.6 Hi ^ gg 

The two-gluon decay mode is handled in a way very similar to Hi 77. The width is [14] 



r{Hi ^ gg) = 1^ (i^f |S^(M^,J|2 + K|, \P^{MH,)f) , 



(29) 
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with the form factors 



S^{MhJ = Cl^q^Fsir,), P^{Mh.) = J] C£,,-^F,(r,) (30) 

q '"'I q "'''i 

and the same expressions given by eqs. ()27p . (j28p for the loop functions. The kinematic variabe 
Tq = iJiff /Amg. In the expression for T{Hi — > gg), the leading-order QCD corrections are known in 
the heavy quark limit {niq ^ M//- ) . These are given by [14] 



95 


7 


T 




97 


7 


T 





(31) 



where is the number of quark flavours for which ruq < rriH^. As a word of caution, it should be said 
that the complete NLO corrections may deviate from what is given here by order 10% or more when 
moving away from the parameter space region where the heavy quark limit is applicable, e.g. at high 
tan P when A couples preferentially to b quarks. 



2.5 Example 1: Tree-level MSSM 

As is well known, the supersymmetric version of the 2HDM requires the Yukawa couplings of the model 
to be of type II. Opposite hypercharges for the Higgs doublets, acting as Z2 quantum numbers, are 
introduced by the requirement that the superpotential is analytic. In addition, most of the parameters 
in the potential are fixed by supersymmetry to values given by the electromagnetic and weak gauge 
couplings. As a result the Higgs sector of the MSSM can be completely described, at tree-level, by the 
two inputs ruA and tan f3. For completeness we give here the potential parameters of eq. ([1]) expressed 
in the electroweak gauge couplings g and g' [15]: 

, , 9^ + 9'^ . 9^ - 9'^ . 5^ 

Ai - A2 - , A3- , A4--Y, 

The tree-level MSSM is unfortunately not a very interesting model since the mass of the lightest 
CP-even Higgs boson always fulfills rrih < mz- Such a light h is in contradiction with experiment [16], 
except under special circumstances when the ZZh coupling - proportional to S/3_a - is small. To study 
the MSSM Higgs sector seriously, it is therefore highly recommended to use one of the specialized tools 
available [17]. These include the higher order corrections to the Higgs masses which are necessary to 
get a phenomenologically viable model. In the MSSM at loop level, X5-7 will in general be non-zero, 
and CV is not necessarily conserved. 



2.6 Example 2: Dark scalar doublet 

A special case of the 2HDM is when the potential has an exact Z2 symmetry in the Higgs basis, i.e. in 
the basis where only one of the two doublets develops a vev [18, 19]. In other words, Ae = A7 = 0, 
771^2 = 0, and tan (3 = oH Such a symmetry results in a conserved parity; we call this ff-parity. 

^Note that this special model cannot be obtained as a continuous tan /3 ^ limit from the general potential. 
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As a consequence of the conserved quantum number, each scalar mass eigenstate originates from a 
particular doublet, and no mixing takes place between the doublets. Assigning a positive if-parity to 
the right-handed fermions, couplings to the scalar doublet without a vev (given H = —1) are forbidden. 
The doublet which does not couple to fermions is termed the dark scalar doublet, or sometimes the 
inert Higgs doublet. Note however that the scalars originating from this doublet still couples to the 
EW gauge bosons. 

The dark scalar doublet model has some interesting phenomenological consequences: First of 
all, particles with negative //-parity have to be produced in pairs, which voids some of the Higgs 
production mechanisms at the LHC. Second, since the lightest H-odd particle is stable, it provides an 
excellent candidate for scalar dark matter [19]. It has also been argued that this model allows for a 
heavier "SM-like" Higgs boson compatible with the EW precision data [18]. 



3 Constraints 

3.1 Positivity of the potential 

For the vacuum configuration to be stable the Higgs potential, given by eq. must be positive 
in all field space directions for asymptotically large values of the fields [20]. At large field values, 
the potential is dominated by the quartic terms. These can be conveniently parametrized [21] by 
introducing 

|^>i|=rcos7, |<i>2| = r sin7, , = pe"^, (32) 

l^il |^2| 

with 7 G [0, vr/2], p £ [0, 1], and 6 G [0, 27r). Dropping the common factor r'*, the quartic terms of the 
potential thus become 

V4=-Aicos 7H — A2 sin 7 + A3 cos 7 sin 7 + A4P cos 7 sin 7 

2 2 (33) 

+ Asp^ cos^ 7 sin^ 7cos(20) + [Ag cos^ 7 + A7 sin^ 7] 2/> cos 7 sin 7 cos 9. 

Minimizing the expression in eq. ()33p , a number of necessary conditions for V4 to be positive everywhere 
can be derived: 

Ai > 0, A2 > 0, A3 > -^/X^. (34) 
When Ag = A7 = the additional condition 

A3 + A4 - IA5I > (35) 

together with those of eq. (I34p , provide necessary and sufficient conditions for positivity. When either 
As / or A7 7^ the modified condition 

A3 + A4-A5 > -TaIA^ (36) 

applies instead. In addition to the conditions given by eqs. ([M|) and ([36|) . three more complicated 
expressions arising from the boundaries cosO = ±1 and p = 1 are needed in this case. For brevity we 
have omitted these expressions here, but they are included in the program. 
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3.2 Tree- level unitarity 



In addition to the constraints on the parameters of the potential obtained from positivity, one can also 
obtain limits by requiring tree-level unitarity for the scattering of Higgs bosons and the longitudinal 
parts of the EW gauge bosons [22] . 

As is well known, the complete all-order scattering (S) matrix has to be unitary. Essentially there 
are two different ways of achieving this. In a weakly coupled theory higher order contributions to the 
S'-matrix typically become smaller and smaller. If the theory on the other hand is strongly coupled, 
the individual contributions may be arbitrarily large. It is then only in the sum that they cancel 
and unitarity is respected. For a weakly coupled theory it is natural to require that the 5-matrix is 
unitary already at tree-level. Necessary and sufficient conditions on the eigenvalues of the 5-matrices 
to achieve this in the general 2HDM have been given in [23] . Expressed in the parameters Aj of the 
general potential, the matrices S(y,t,) for scattering of states with specific total hypercharge Y and 
weak isospin E, are given by 

Ai As V^Xq 
As A2 V2A7 
v^Ae v^Ay A3 -I- A4 



167r5, 



(2,1) 



(2,0) 



167r5, 



(0,1) 



A3 — A4 



Ai 


Ai 


A6 


A4 


A2 


At 


Ae 


A7 


A3 


Ae 


A7 


As 



As 
At 
As 
A3 / 



(37) 



167r5, 



(0,0) 



3Ai 
2A3 + A4 
3A6 

3A6 



2A3 + A4 
3A2 
3A7 

3A7 



3A6 
3A7 
A3 + 2A4 
3A5 



3A6 
3A7 

3As 
A3 + 2A4 



To respect tree-level unitarity, the eigenvalues Lj of these matrices should fulfill |Lj| < IGtt, which 
corresponds to saturating the ^-matrix with the tree-level contribution. Alternatively one can make 
a partial-wave expansion of the scattering amplitudes for the different channels as discussed in [1] and 
put limits on the partial wave amplitudes such as |Re(aj)| < 1/2. The requirement |Re(ao)| < 1/2 for 
the J = amplitude then corresponds to |Lj| < Sir. Finally, one can also impose (possibly harder) 
constraints on the parameters of the potential based on arguments of perturbativity, by demanding 
that the quartic Higgs couplings fulfill {CHiHjH^Hil < 47r or similar constraints. (In principle one 
could also imagine a mixed situation where the interaction in some channels is strong and in others 
perturbative [23].) 



3.3 Oblique parameters 

The electrowcak oblique parameters S, T, U [24] constitute a sensitive probe of new physics coupling to 
the EW gauge bosons. Values on these parameters are constrained by the precision measurements at 
LEP (see [25]). In general, the contribution to the oblique parameters from the extra SU{2) doublet in 
the 2HDM is small, since scalar doublets (or singlets) do not break the custodial symmetry [26] which 
protects the tree-level relation p = Mw /{Mz cos 6w) = 1- However, large mass splittings among the 
new Higgs states can still induce sizeable contributions at the loop level. 
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Prom the general result of [27], we determine the 2HDM contribution to the oblique parameters. 
For brevity, only the expressions for S, T and U are given here, but the full set of six oblique 
parameters [28], i.e. also including V, W and X, are available in 2HDMC. Starting with T, this parameter 
is related to the commonly used parameter po = P/PSM through po = 1/(1 — aT). It thus encodes the 
departure from the SM value oi po = 1. The 2HDM contribution to T is given by 
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647r2 



m 



w I k=l 



k=l 



+ 3^|C3-fc|' \f {ml,mlJ-F{mlr,mlJ 



k=i 



(38) 



F {ml, ml^ref) - F (m^, m^^,gf) 



where the masses of the neutral 2HDM Higgs bosons arc denoted mn,, {k = 1, 2, 3), and nih^^ci is the 
mass of the SM Higgs boson. The contribution of the latter is subtracted in order to isolate the new 
physics contribution. The function F is given by 



, x + y xy , X 

F{x,y) = ^ ^In-, 

2 x-y y 



(39) 



and the coupling structure of the EW gauge bosons to a pair of Higgs bosons has been summarized 
in the vector 

(7= |c;3-a,S/3_a,l}- (40) 

Before giving the expressions for the S and U parameters, two additional auxiliary functions, G and 
G, are needed. The first of these is defined as 



x'^ + — (x — y) 



x-y Q 

where t = x + y — Q,r = Q'^ — 2Q{x + y) + {x — y)^, and 

t- 



+ 



3Q2 



+ (41) 



f{t,r) = { 







t + y/r 



2\/^ arctan 



r > 0, 
r = 0, 
r < . 



(42) 



The second function G is given by 

G{x, Q) = G{x, Q, Q) + 12G(x, Q, Q), 

where 



G{x,y,Q) = -2 + 



X — y x + y 
Q x-y^ 



ln^ + i/(.,.) 



(43) 
(44) 
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Having defined all these functions, the S parameter can then be written simply as 



S = 



3847r2c2^ 



k=l 



i^w -CwfG {m'jj+ , m'jj+ , m|) + ^ | Q |^ G {m%^ , m\, m|) 

k-- 

3 

- 2 In + ^ In - In m^^^^f 



(45) 



fc=i 



k=l 



where sw = sin^^, cw = cos 6w, and Ow is the weak mixing angle. Finally, the U parameter is given 

by 



U 



2 f 3 

^\Gkf G {mjj+,mjj^,mlr) - [s^, - c^]^ G {m]j+,m]j+,m%) 
I k=l 

2 

k=l 
2 

+ Yl G (m^^, m^) - G {m%^,m\) - G (m^ m^) + G (m^^.^f, m|) 



ifc=i ; 

(46) 

The current experimental limits on 5, T and C/ can be found in PDG [25]. It should also be noted 
that typically the S and T parameters cannot be fitted separately from the mass of the SM Higgs 

boson, nih^j-cf. Therefore, when comparing with experimental limits, care has to be taken to use the 
same value of TO/i,ref when calculating the oblique parameters in the 2HDM at hand, as was done in 
the SM fit one is comparing to. 



3.4 Anomalous magnetic moment of the muon 

The extended Higgs sector of the 2HDM contributes to the anomalous magnetic moment of the muon, 
(^^i = {g^i ~ 2)/2. Since the experimentally measured value for shows a deviation of So" from the 
SM value [29], this provides an interesting test of new physics models. An overview of the current 
status, with updated values on from both experiment and theory, can be found in [30]. 

At one loop, both the charged and neutral Higgs bosons of the 2HDM add to the vertex correction, 
but these corrections are generally small compared to the observed deviation [31]. More interesting are 
the dominant positive contribution at two loops, which is of Barr-Zee type [32] from diagrams involving 
a pseudoscalar Higgs, a photon propagator, and a heavy fermion in the upper loop. The large fermion 
mass compensates with a factor m j / for the loop suppression. Contributions of similar type, albeit 
supressed by orders of magnitude, exist also from combinations of HiZ and H^W^. From [33], we 
get the 2HDM contribution Sa^ to from a photon Barr-Zee diagram with a neutral Higgs and a 
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heavy fermion / in the loop 



^<=i^'v.<3j!!;;|£2r3VpN(w) 

where y/ = rnj/m^^_, Nc = 3(1) for quarks (leptons), Qf the electric charge, and t/ the third 
component of the weak isospin (T3 = —1/2 for a down-type fermion). With slight abuse of notation, 
we have used K^^p^, etc. to denote the appropriate diagonal elements of the corresponding Yukawa 
matrices indexed with U, D, L above. Finally, the functions / and g are given by the integrals 

,(.)=ra,l_Mlz^,,£(l^, ,48) 

Jo x{l-x)- z z 



and 



g{z) = I ' dx--^ In ^i^. (49) 

x{l — x) — z z 



In 2HDMC, the dominant two-loop contributions according to eq. ()47p have been included for all 
fermions. For completeness, the sub-dominant contribution from the one-loop corrections has also 
been included. 



3.5 Higgs searches at colliders 

Since no Higgs boson was observed so far either at LEP-II or the Tevatron, stringent limits have 
been derived on the masses and couplings allowed in the Higgs sector, both within and beyond the 
SM. Rather than to implement a full analysis of the collider constraints in the general 2HDM, we 
provide interfaces for 2HDMC to use the existing codes HiggsBounds [34] or NMSSMTools [35] to check 
consistency with the existing limitsU HiggsBounds includes most search channels for neutral Higgs 
bosons at LEP, and in addition a large number of channels explored by the Tevatron experiments. 
NMSSMTools contains most of the relevant search channels for neutral Higgs bosons at LEP. For an 
up-to-date list of the channels included, as well as a description of the statistical methods employed 
to combine the results from several channels, we refer to the original publications [34, 35] describing 
these codes. 

For the charged Higgs boson, we have implemented the DELPHI limits from the general 2HDM 
search [36]. The two limits which could be used are both from the pair production process e^e^ 
H~^H~, followed by decays in the — > t^Ut and cs channels. These two modes are 

particularly important in the type II 2HDM. 

''Consult section [5] of this manual for more information on how to install 2HDMC with HiggsBounds or NMSSMTools 
support. 
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Parameter 


Value 


Parameter 


Value 




127.91 


Gp 


1.16637 X 10-5 GeV-2 


Osimz) 


0.1176 


niz 


91.1876 GeV 



Table 2: Default values used in 2HDMC for the Standard Model gauge parameters. 



Parameter 


Value 


Parameter 


Value 


Parameter 


Value 


ruu 





rrid 







0.511 X 10-^ 


rric 


1.77 




0.1 




0.10566 


mt 


171.2 


rUb 


4.8 


m-j- 


1.777 



Table 3: Default values used in 2HDMC for the fermion pole masses. All values are in GeV. 

4 2HDMC code structure 

4.1 General description 

This section of the manual describes the structure of the 2HDMC code in terms closely related to the 
physics. A more technical documentation of the 2HDMC classes and functions, such as lists of function 
arguments, are given in the full class documentation distributed with the code0 

2HDMC is built as a modular C++ code along standard principles for object-oriented programming. 
Every program should therefore contain one (or more) THDM object(s) representing the 2HDM(s) under 
study. Since the 2HDM is an add-on to the SM, each THDM contains its own instance of the SM class. 
This allows for flexible specification of the underlying SM parameters. 

The THDM class itself provides member functions to operate on the 2HDM in various ways, such 
as changing the basis or the potential parametrization. It further contains the specification of the 
Yukawa sector, and functions to return the couplings of the 2HDM Higgs bosons. This class also 
handles input/output to file. To calculate the decay widths and branching ratios of the 2HDM Higgs 
bosons, the class DecayTable is provided. Objects created from this class will operate on a THDM. The 
same is true for the Constraints class, which contains functions to calculate the various observables 
discussed in section [3l 

4.2 Physics input 

The Standard Model parameters in 2HDMC are contained in objects of the SM class. Default values for 
all parameters (taken from PDG [25]) are stored as static members in this class. The default values 
are tabulated in table [2] (for the gauge sector parameters) and table [3] (fermion masses). To obtain the 
SM object which contains the SM parameters for an instance of THDM, the function get_SM should be 
used. The parameters can then be changed by invoking different functions, e.g. to set the top quark 
mass set_qmass can be used. For a complete list we refer to the SM class documentation. Finally, 
to use the new parameters, the function set_SM is called for the THDM with the modified SM object as 
argument. An example of how to set mt is given in the Demo program. 

In the EW sector, the SM is specified by three parameters: a, Gp, and mz, from which mw and 
the weak mixing angle sin^ 9w = 1 — "ifv/'^z determined. From these basic parameters, values 

' |http : //www 7isv .'uu'.'se/thep/MC/2HDMC | 
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THDM function Required parameters 



Comment 



set_param_gen Ai, A2, A3, A4, A5, Ag, A7, mfg, tan/3 General potential, eq. ([T|) 

set_parainJiiggs Ai, A2, A3, A4, A5, Ag, A7, mH± Higgs basis, tan/3 = 

set_parain_hhg Ai, A2, A3, A4, A5, Ag, tan/3 Z2-symmetry only softly broken 

set_param_phys rrih, mn, rriA, mu±, -S/^-q, Ag, A7, m^2' tan/3 Physical masses 

Table 4: Input parameters for the Higgs potential in different bases and parametrizations. 



for the supplementary EW quantities are computed according to 

v = [y2Gp) ' , g = , g = tant^iy, e = smfc'iy. (50) 

V V V 

The 2HDMC user can enter the Higgs potential parameters in a number of different parametrizations 
using member functions of THDM. A full list of the available parametrizations, and their corresponding 
parameters, is given in tableUl For each function set_parain_X listed in this table, there is a correspond- 
ing function get_param_X which can be used to retrieve the current model in that parametrization. 

A few additional comments about table H] are in place. The number of required input parameters 
depends on the assumptions underlying the specific parametrization. For instance, in the Higgs basis 
tan /3 = and need not be specified. Another example is given by the potential in the Higgs Hunter's 
Guide [1] which was originally introduced in [37]. This parametrization only allows for soft breaking 
of the Z2 symmetry, equivalent to forcing Ag = A7 = in the general potential. Note the different 
meaning of Aj in this parametrization compared to the general case (the relation between the two can 
be found, for example, in Appendix A of [2]). In the physical basis, the Higgs masses are given as input 
together with the invariant S/3_a- The remaining parameters in this parametrization are of course not 
more physical than in any general potential. Since Ag,A7, and play a role in the Higgs self- 
interactions, their effects appear for instance in the loop-mediated decay — > 77. In the conventions 
used, the allowed range for a is fixed by the requirements < /3 < 7r/2, and —tt/2 < f3 — a < Tr/2. 

Having specified the Higgs potential, we can move on to specify the Yukawa sector in order to have 
a complete model. This is only required when Higgs couplings to fermions are going to be considered. 
Several THDM functions, listed in table O exist for specifying the Yukawa sector according to eq. (jlip . 
Non-diagonal elements in the p matrices are allowed, but will lead to FCNCs. There are presently 
no checks against limits on FCNC included in the program. It is important to note that k is defined 
internally using pole masses. When related to k, the elements of p should therefore also be specified 
using pole masses. To resum large logarithms, the flavour-diagonal couplings Hiqq are rescaled with 
rnq{mHi)/mq°^'^ . Specifying p matrices directly related to e.g. the MS masses can therefore lead to 
unwanted consequences, such as a broken Z2 symmetry in the Yukawa sector. 

Two special cases of the 2HDM have been implemented as special input functions in THDM. These 
are described in table [H What distinguishes these are that they set both the potential parameters 
and the Yukawa sector using a single function. Not surprisingly, set_MSSM sets the 2HDM parameters 
according to the tree- level MSSM relations (see section r2.5p . with the Yukawa sector of the MSSM 
(type II). Again we emphasize that the tree-level MSSM in only included for reference, and that it 
is not intended to be used in any physics analysis. The other special function is set_inert, which 
sets the parameters of the dark scalar doublet model (section 12. 6p . The convention chosen for the 
parameters of this model is that of [38]. Only one doublet is given Yukawa couplings in this case. 
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THDM function 


Arguments 


Comment 


set_yukawas_type 

set_yukawas_up 

set_yukawas_up 


Type (1-4) 

V U U 
ruu^ rcc^ rtt 


Universal Yukawa type (see table [1]) 
Elements of diagonal matrix 
Elements of full p^ matrix 



Table 5: Functions which are used to specify the up-type quark Yukawa sector. Similar methods 
exist for the down-type quark and lepton Yukawas. Note that set_yukawas_type specifies all sectors 
simultaneously. Yukawa matrices related to quark masses should always be specified in terms of the 
quark pole masses. 



THDM function 


Required parameters 


Comment 


set_MSSM 
set_inert 


ruA, tan/3 

mf^, m^, m^, A2, A3 


Tree-level MSSM potential, Yukawas type II 
Dark scalar (DS) doublet has no vev, no Yukawas 



Table 6: Special functions used for giving input parameters to the Higgs potential, including specifi- 
cation of the Yukawa sector. 



4.3 Program output 

The main output which is possible to get from 2HDMC can be summarized as follows: 

• Potential parameters in different parametrizations and bases 

• Higgs masses and couplings 

• Higgs decay widths and branching ratios 

• Properties of interest for constraining the 2HDM parameter space 

It is very simple to use 2HDMC for converting between different parametrizations of the Higgs potential. 
Simply give the parameters using one of the set_param functions, followed by a get_param to get 
another representation of the same model. The basis can also be changed in this way. Invoking 
recalc_taii_beta, the model parameters are recalculated for a new value of tan/?. Of course this 
does only work for the {7(2)-invariant Higgs potential, not for the Yukawa sector. To write out the 
parameters in a certain model to stdout, a number of THDM functions starting with print_param are 
available. 

When the Higgs masses are not used as input, they are calculated from the tree- level Higgs po- 
tential. The mass values can be accessed either through get_paraiii_phys, or using getJimass which 
takes an index for the Higgs boson as argument. For every method where this is the case, we use the 
enumeration {h = 1, H = 2, A = 3,H^ = 4}. The couplings involving one or more Higgs bosons are 
available through the functions listed in table [71 

To calculate the partial decay widths of the 2HDM Higgs bosons into various channels, the class 
DecayTable is used. The functions available through this class are collected in table [HI To toggle 
QCD corrections and loop-induced decay modes {Hi gg,J'j) off/on, the method set_QCD can be 
used. 

As discussed in section [Sj there exist theoretical constraints on the 2HDM parameters. Functions 
to check that these are satisfied are provided in the Constraints class. This class also provides the 
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THDM function 



Coupling 



Fcynman rule Convention 



ge t _c oup 1 ingJidd 


Hi djdk 


Cs + 






get_couplingJiuu 


Hi UjUk 


Cs + 


IbCp 




get_couplingJill 


Hi 


Cs + 


IbCp 




get _coupl ingJiud 


H~^djUk 


Cs + 


75 Cp 




get_coupling_hln 




Cs + 


75 Cp 




get_coupling_vvh 


HiVjVk 






C in GeV 


get_coupling_vhh 


HiHjVk 


C{pi- 


-pjr 


All momenta ingoing 


get-couplingjihh 


HiHjHk 


C 




C in GeV 


get_coupling_vvhh 


HiHjVkVi 


Cgt"" 




Single W'^ momenta ingoinj 


get_couplingJihhh 


HiHjHkHi 


C 







Table 7: Functions which arc used to calculate the Higgs boson couplings (denoted C in the table) . For 
a purely neutral current interaction, i = 1 . . . 3. When both neutral and charged currents are possible, 
i = 1 . . . 4 with as i = 4. Vector indices go from 1 to 3, corresponding to 7, and W'^ . Fermion 
family indices go from 1 to 3. The flavour-diagonal couplings to quarks assume renormalization of 
quark MS masses at the Higgs mass scale. 



DecayTable function Width Comments 



get .gammaJidd 


V{Hi^ 


djdk) 


NLO QCD corr. 


get _gamma Jiuu 


T{Hi^ 


UjUk) 


NLO QCD corr. 


get _gainmaJill 


T{Hi^ 


ijik) 




get -gammaJiud 


T{H+ - 


^ djUk) 


LO QCD corr. 


get .gamma Jiln 


T{H+ - 


- i^i^k) 




get -gamma Jivv 


r{Hi^ 


V^*Wj) 




get _gainma_hvh 




V}*^H,) 


Widths for W+R- and W H+ added 


get _gamma_hhh 


V{Hi^ 


HjHk) 




get -gamma Jigaga 




77) 


Effective operator 


get _gammaJigg 




99) 


Effective operator, LO QCD corr. 



Table 8: Functions which calculate the decay widths of Higgs boson Hi, i = 1 . . .4. Fermion family 
indices go from 1 to 3, whereas gauge boson indices go from 2 to 3, meaning in turn Z and W^. The 
decay widths for H and H V^*''H also include the decay through a virtual gauge boson. 
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Constraints function Description 



Returns 



check_positivity 

check_unitarity 

check_perturbativity 

check_masses 

oblique_param 

delta_amu 



Checks potential positivity true/false 

Checks tree- level unitarity true/false 

Checks perturbativity true/false 

Checks the collider constraints true/false 

Calculates EW oblique parameters {S, T, U, V, W, X) values 

Calculates the 2HDM contribution to o„ value 



Table 9: Functions which are available in the Constraints class. Functions with names beginning with 
check_ give boolean true for models that comply with the required condition, and false otherwise. 
The other functions return the 2HDM contributions for comparison to experimental limits. 



Generic basis 



MINPAR 






3 


tan/3 


Set to zero for Higgs basis 


11-17 


A1-A7 




18 


mfa 


Generic basis only 


MASS 






37 


mH+ 


Higgs basis only 



Table 10: LesHouches input parameters for the generic and Higgs bases 



calculation of the 2HDM contribution to the other observables discussed in section [3l Table [9] lists the 
various functions that can be of interest when using 2HDMC to constrain the 2HDM parameters. 

Another set of 2HDM constraints is given by various flavour-physics observables, such as rare 
decays of B-mesons. No flavour observables are included in 2HDMC, but instead the code can be easily 
interfaced with Superlso [39] . 

4.4 LesHouches Interface 

The input to 2HDMC can also be specified through a file interface inspired by the LesHouches format [40]. 
To distinguish it from parameter files for the MSSM, the file must contain 

Block MDDSEL 
10 

to select the 2HDM. The Higgs potential input and basis choice can be given in three formats, corre- 
sponding to the general basis, the Higgs basis, and physical basis as described above. For the general 
and Higgs bases, the relevant blocks and parameter numbers that need to be specified are given in 
table [TOl while those for the physical basis are given in table [TTJ 

When tan/3 7^ 0, the type of Yukawa couplings can also be specified using MINPAR(24) to select 
between types 1-4. The LesHouches interface cannot be used to specify more generic Yukawa couplings 
at this point. A sample input file named LH_input is included in the 2HDMC distribution. This file can 
be run with the example program CalcLH. 2HDMC can also produce LesHouches- formatted output by 
calling the THDM function writeXesHouches with the filename specified. For an example of how such 
an output file might look, see appendix iBl 
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Physical basis 



MINPAR 
3 

16-17 



tan /3 



20 



18 




MASS 



25 
35 
36 
37 



ruH 



ruA 



Lightest CP-even Higgs 
Heaviest CP-even Higgs 
CP-odd Higgs 
Charged Higgs 



Table 11: LesHouches input parameters for the physical mass basis. 



4.5 MadGraph/MadEvent model 

A USRMOD model for MadGraph/MadEvent [41] is provided with 2HDMC to facilitate event generation in 
the 2HDM. In the MGME subdirectory of the 2HDMC distribution, a directory called 2HDMC contains the 
model. This directory should be placed in the Models directory of MadGraph/MadEvent. Using the 
standard output function write_LesHouches with the couplings argument produces model output in 
the format of a param.card. All couplings listed in tabled except the four-Higgs vertices, are included 
in the 2HDMC model. At this point the model cannot handle flavour-changing neutral currents. To use 
tree- level widths consistently, the 2HDMC QCD corrections to Higgs decay widths into fermions and the 
decays Hi — > gg,^^ can be swtiched off using the set_QCD function of DecayTable. 

One advantage with our model over the "default" 2HDM model in MadGraph/MadEvent [41] is the 
treatment of couplings and decay widths. All input to our model is computed with 2HDMC to ensure 
internal consistency. The Hiqq couplings (evaluated at the correct Higgs mass scale) are then input 
to the model. With the default model, these couplings are computed in the model itself, whereas the 
decay widths are computed externally by a calculator. An obvious advantage with the existing model 
is that it allows for simulation of processes that violate CP, which our model does not. To use the 
output of 2HDMC together with the default 2HDM model should require only minor modifications of 
the LesHouches file format used. 

5 Downloading and compiling the package 

The latest version of the 2HDMC code can be obtained from 

http : //www . isv . uu . se/thep/MC/2HDMC . 

In addition, this manual and the complete class documentation, including the specification of the 
required input and output of various functions, can be downloaded from the same URL. To compile 
2HDMC, the following tools are required: 

• A C++ compiler, such as GCC 

• Common build utilities: make etc. 
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• The GNU Scientific Library (GSL); available from |http : / /www . gnu . org/sof tware/gsl/ 

The 2HDMC library has been successfully compiled using GCC 4 under Linux. However, there should 
be no difficulties in using GCC 3 or another compiler and/or architecture. Please report to the authors 
your experiences from compiling and running 2HDMC on another platform. Unpacking the main tarball, 
a directory called 2hdmc is created. To start the compilation, enter this directory and run 

$ make 

to compile the library and some example programs (file names starting with Calc) that demonstrate 
the use of 2HDMC in standalone applications. The installation can be verified by running the Demo 
program and comparing the output to appendix [Aj To compile the 2HDMC library only, run 

$ make lib 

while a main program (Program, cpp, located in 2hdmc/src) is compiled with the 2HDMC library using 
the command 

$ make Program 

We encourage the user to examine the Makefile for details. 

As described in section 13.51 there are two possibilities to access constraints from collider experi- 
ments in 2HDMC. These can be used either independently or in combination, and neither is required to 
run 2HDMC. However, if you are using any of these codes with your 2HDMC project, make sure to cite 
the relevant publication. 

HiggsBound^ [34] provides limits from both LEP and the Tevatron experiments. To use HiggsBounds, 
the library libHB . a should be placed in the 2hdmc/lib directory, and in addition two lines should be 
uncommented in the 2HDMC Makefile. The interface has been tested to run with HiggsBounds 1.0.3. 

NMSSMTool^^ [35] includes limts from LEP only. To use NMSSMTools, the library nmhdecay.a 
and the experimental data provided in the distribution is needed. The library should be put in the 
2hdmc/lib directory, and the experimental data must be accessible in . ./EXPCON when running a 
program compiled to use NMSSMTools. To compile 2HDMC with NMSSMTools support, three lines should 
be uncommented in the Makefile. Version 2.2.0 has been tested to work. 
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A Sample output 



2HDMC ships with numerous example main programs. Compihng and running Demo . cpp should produce 
the following output to stdout: 



$ ./Demo 



* * 

* 2HDMC - Two-Higgs-Doublet Model Calculator * 

* http://www.isv.uu.se/thep/MC/2HDMC * 

* Version 1.0.3 * 

* Compiled on Jun 22 2009 * 

* * 



2HDM parameters in physical mass basis: 



m_h 


80 


00000 


m_H 


200 


00000 


m_A 


140 


00000 


m_H+ 


160 


00000 


sin(b-a) 





20000 


lambda_6 





00000 


lambda_7 





00000 


ml2'2 


800 


00000 


tan (beta) 


10 


00000 



2HDM paramet 
lambda_l : 
lambda_2 : 
lambda_3 : 
lambda_4 : 
lambda_5 : 
lambda_6 : 
lambda_7 : 
ml2~2: 

tan (beta) : 



s in generic basis: 

2.25553 

0.61585 

2.29636 
-0.38796 
-0.19002 

0.00000 

0.00000 
800.00000 
10.00000 



2HDM paramet 
Lambda. 1 : 
Lambda_2 : 
Lambda_3 : 
Lambda_4 : 
Lambda_5 : 
Lambda_6 : 
Lambda_7 : 
m_Hp: 



s in Higgs basis: 
0.63763 
2 . 24484 
2.29082 
-0.39350 
-0.19556 
-0.10861 
-0.05374 
160.00000 
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Constraints : 

Potential stability: OK 

Tree-level unitarity: OK 

Perturbativity : OK 

Mass constraints: OK 



Oblique S 


6 


77200e- 


-02 


Oblique T 


-1 


62816e- 


-02 


Oblique U 


2 


95734e- 


-04 


Oblique V 


7 


70150e- 


-03 


Oblique W 


1 


02287e- 


-03 


Oblique X 


7 


70645e- 


-04 


Delta_rho 


-1 


27289e- 


-04 


Delta_amu 


3 


78746e- 


-11 



Total width for h : 1.707e-01 GeV 

Total width for H : 5.212e-(-00 GeV 

Total width for A : 2.979e-01 GeV 

Total width for H-^: 4.083e-02 GeV 

B LesHouches output 

The same Demo program discussed in appendix [A] produces the LesHouches file Demo_out . lha with 
the following content: 



################################################################## 



# # 

# Two-Higgs Doublet Model LesHouches-card # 

# Produced by 2HDMC: www.isv.uu.se/thep/MC/2hdmc # 

# # 



################################################################## 
Block MODSEL # Select Model 
10 # 10 = THDM 

Block SMINPUTS # Standard Model inputs 



1 1.27910000e+02 # l/alpha_em(MZ) SM MSbar 

2 1.16637000e-05 # G Fermi 

3 1.17600000e-01 # alpha_s(MZ) SM MSbar 

4 9.11876000e+01 # MZ 

5 4.24693805e-^00 # mb(mb) 

6 1.75000000e+02 # mt (pole) 

7 1.77700000e+00 # mtau(pole) 



Block MINPAR # Model parameters 

# Parameters for general potential in generic basis 
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3 l.OOOOOOOOe+01 # tan (beta) 

11 2.25553058e+00 # lambda. 1 

12 6.158525366-01 # lambda_2 

13 2.296360456+00 # lambda_3 

14 -3.879615238-01 # lambda_4 

15 -1.900219716-01 # lambda_5 

16 0.000000006+00 # lainbda_6 

17 0.000000006+00 # lambda_7 

18 8.000000006+02 # m_12-2 

20 2.000000006-01 # sin (beta-alpha) 

21 9.797958976-01 # cos (beta-alpha) 
24 2 # Yukawas Type 

Block MGCKM # CKM elements 

1 1 9.741900006-01 # Vud 

1 2 2.257000006-01 # Vus 

1 3 3.590000006-03 # Vub 

2 1 2.256000006-01 # Vcd 
2 2 9.733400006-01 # Vcs 

2 3 4.150000006-02 # Vcb 

3 1 8.740000006-03 # Vtd 
3 2 4.070000006-02 # Vts 
3 3 9.991330006-01 # Vtb 

Block MASS # Mass spectrum (kinematic masses) 

# PDG Mass 

1 O.OOOOOOOOe+OO # Md 

2 O.OOOOOOOOe+OO # Mu 

3 l.OOOOOOOOe-01 # Ms 

4 1.770000006+00 # Mc 

5 4.800000006+00 # Mb 

6 1.750000006+02 # Mt 
11 5.109989186-04 # Me 
13 1.056583676-01 # Mmu 
15 1.777000006+00 # Mta 

23 9.118760006+01 # MZ 

24 7.982572886+01 # MW 

25 8.000000006+01 # Mhl, lightest CP-even Higgs 

35 2.000000006+02 # Mh2, heaviest CP-even Higgs 

36 1.400000006+02 # Mh3, CP-odd Higgs 

37 1.600000006+02 # Mhc 

Block alpha # Effective Higgs mixing parameter 
1.269769756+00 # alpha 

# PDG Width 

DECAY 6 1.556882946+00 # top decays 

# BR NDA IDl ID2 
9.988810186-01 2 24 5 
2.258826326-08 2 37 1 
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DECAY 25 
# 



DECAY 35 
# 



DECAY 36 
# 



4.91600304e-07 


2 


37 


3 


1.11846795e-03 


2 


37 


5 


25 1. 706694886 


1-01 


# hi 


decays 


BR 


NDA 


IDl 


ID2 


1.06791622e-03 


2 


3 


-3 


6.75570041e-05 


2 


4 


-4 


9.06817540e-01 


2 


5 


-5 


7.40019764e-09 


2 


11 


-11 


3.16378237e-04 


2 


13 


-13 


8.92259685e-02 


2 


15 


-15 


3.28604940e-06 


2 


22 


22 


2.50134614e-03 


2 


21 


21 


35 5.21246842e+00 


# h2 


decays 


BR 


NDA 


IDl 


ID2 


8.17946511e-06 


2 


3 


-3 


4.80480461e-05 


2 


4 


-4 


6.087365656-03 


2 


5 


-5 


5.838631086-11 


2 


11 


-11 


2.49619285e-06 


2 


13 


-13 


7.057326166-04 


2 


15 


-15 


2.391950486-05 


2 


22 


22 


6.817634126-02 


2 


23 


23 


1.94337426e-01 


2 


24 


-24 


1.728290556-04 


2 


21 


21 


7.304333096-01 


2 


25 


25 


2.595945806-06 


2 


23 


36 


8.78581628e-07 


2 


24 


-37 


8.785816286-07 


2 


37 


-24 


36 2.97857573e 


1-01 


# h3 


decays 


BR 


NDA 


IDl 


ID2 


1.14103374e-03 


2 


3 


-3 


6.834551116-06 


2 


4 


-4 


8.992125616-01 


2 


5 


-5 


8.055098006-09 


2 


11 


-11 


3.44380190e-04 


2 


13 


-13 


9.737909656-02 


2 


15 


-15 


3.622405626-07 


2 


22 


22 


1.073459816-03 


2 


21 


21 


8.42263714e-04 


2 


23 


25 



# Charged Higgs decays 





BR 


NDA 


IDl 


ID2 


4, 


.699222146-04 


2 


2 


-3 


9. 


, 17812917e-05 


2 


2 


-5 


2, 


.753482416-06 


2 


4 


-1 


8, 


.788706076-03 


2 


4 


-3 


1, 


.226186906-02 


2 


4 


-5 
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6.71516661e-08 
2.87094082e-03 
8.11866349e-01 
1.63559670e-01 
8.79405283e-05 



2 -11 12 

2 -13 14 

2 -15 16 

2 24 25 

2 24 36 
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